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Abstract 

According to Dirac's theory of the positron, an electromagnetic field 
tends to create pairs of particles which leads to a change of Maxwell's 
equations in the vacuum. These changes are calculated in the special case 
that no real electrons or positrons are present and the field varies little 
over a Compton wavelength. The resulting effective Lagrangian of the 
field reads: 
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The expansion terms in small fields (compared to £ ) describe light-light 
scattering. The simplest term is already known from perturbation theory. 
For large fields, the equations derived here differ strongly from Maxwell's 
equations. Our equations will be compared to those proposed by Born. 
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The fact that electromagnetic radiation can be transformed into matter and vice 
versa leads to fundamentally new features in quantum electrodynamics. One 
of the most important consequences is that, even in the vacuum, the Maxwell 
equation have to be exchanged by more complicated formulas. In general, it 
will be not possible to separate processes in the vacuum from those involving 
matter since electromagnetic fields can create matter if they are strong enough. 
Even if they are not strong enough to create matter they will, due to the virtual 
possibility of creating matter, polarize the vacuum and therefore change the 
Maxwell equations. This polarization of the vacuum to be studied below will 
give rise to a distinction between the vectors S, € on the one hand and 2), 5) 
on the other, where 
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The polarizations 03 and 071 can be arbitrary functions of the field strengths at 
the same place, their derivatives, and the field strengths in the surroundings of 
the observed position. If the field strengths are small (which means, as we shall 
see small compared to e 2 /Sc-times of the field strength at the boundary of the 
electron), then 03 and 971 can be approximately considered as linear functions 
of £ and 03. In this approximation, Uehling 1 and Serber 2 have calculated the 
modifications of Maxwell's theory. Another interesting case is obtained by not 
assuming small field strengths but instead slowly varying fields (i.e., the fields 
are nearly constant over the length h/mc). Then one obtains 03 and 971 as func- 
tions of <B and 03 at the same position. The derivatives of l£ and 03 do not 
appear in the approximation. The expansion of 03 and 971 in powers of (£ and 03 
will contain only odd powers, as will be seen in the calculation. The expansion 
terms of third order are phenomenologically related to light-light scattering and 
are already known 3 . The goal of this work is to find the functions 03 ((£,03) 
and 97l(£, 03) for slowly varying field strength. It is sufficient to calculate to 
energy density U(<£, 03). From the energy density one can derive the fields using 
the Hamiltonian method; one introduces the Lagrangian £((£, 03) and obtains 

Lagrangian is determined by (3) and 27 and fj are determined by (2). Due to 
relativistic invariant the Lagrangian can only depend on the two invariants 
£ 2 - 03 2 and (<*B) 2 (compare EJ). The calculation of U{<£, 03) can be reduced 
to the question of how much energy density is associated with the matter fields 
in a background of constant fields £ and 03. Before solving this problem, thc 
mathematical scheme of the positron theory 4 will be presented. This will also 
correct some errors in earlier formulas. 



1 The Mathematical Schema of the Theory of 
the Positron 

Starting point of the theory is thc Dirac density matrix, which is given in naive 
wave theory by 

(x't'k'\R\x"t"k") = ^(n) ^*(x"t"k")ip n (x't'k') (4) 

occupied states 

1 E. A. Uehling, Phys. Rev. 48, 55, 1935. 

2 R. Serber, Phys. Rev. 48 , 49, 1935. 

3 H. Euler, B. Kockel, Naturwissensch. 23, 246, 1935. 

4 W. Heisenberg, ZS. f. Phys. 90, 209, 1934 
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and in this quantum theory of wave fields by 5 



(x't'k'\R\x"t"k") = tP*(x"t"k") i/j(x't'k') 



(5) 



Apart from this matrix an important role is played by the matrix R$, formally 
defined by 



(x't'k'\R s \x"t"k") = U £(n) - £(n) \r n (x"t"k")^ n {x't'k') 



>ccupicd state 



occupied state 



and in the quantum theory by 



(x't'k'\R s \x"t"k'' 



^{x"t"k") iP{x't'k') - ^(x't'k') ^{x"t"k") 



(6) 



(7) 



The matrix R$ as a function of the differences x' x — x'l = x\ and t' — t" = t 
becomes singular on the light cone. If we set 

ct == x, = - ^ = * ; ft - 5L±-5i, (B) 

furthermore the potentials Ao = —A ; Ai = A 1 , and for the Dirac matrices 
a = —ao = \; a 1 — ai, we derive 



W k' \R s \x"k") = u _ ^ + . % , ^ 1 1 



where 6 
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(The sum of repeated Latin indices runs from 1 to 3, and of Greek indices from 
to 3). The integral is to be calculated along the straight line from P' to P" . 

The density matrix r responsible for the behavior of matter is obtained from 
Rs by the equation 

where S is given by 

p" 



8 = e 



a 



X^ X* 



C 



(12) 



Here So is the matrix Rs in field-free and matter-free space, where a, 6, and C 
are defined 7 by 

5 In the earlier paper [1] the values on the right hand side are interchanged 
6 In the earlier paper ^ the exponent has wrongly a negative sign. 

7 In the earlier paper [1] there is a mistake that leads to another value of C. There, the 
character 7 stands for the logarithm of the Euler constant, in contrast to the common usage. 
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where 7 is the Euler constant 7 = 1.781... 

The four-vector for the current density and the energy momentum tensor 
are obtained from r by 



*' ft" 



* S«, t .( f + f,M S -f,r). 

In the quantum theory of wave fields it is useful to expand the wave function 
into an orthogonal system: 
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The operators a n can be represented in the form 

a* - N n A n V n ; a n = V n d n N n > 

where A„ changes the number N n to 1 — N n and V n = Y\ t<n (l — 2N t ). 
We further set 

a n = a* = — V n A K K; a'* = —N' n A«y n ; K = l~N n . 

The Hamiltonian of the whole system reads in these variables: 
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According to the powers of the elementary change e, we distinguish the terms: 
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2 Calculation of the Energy Density in Intuitive 
Wave Theory 

Since the Lagrangian associated with the corrected Maxwell equations is a func- 
tion of the two invariants € 2 — 03 2 and (€B) 2 , it is sufficient to calculate the 
energy density of the matter field as a function of two independent field quan- 
tities. For example, it will be sufficient to calculate the energy density in a 
constant electric and a constant parallel magnetic field. In these constant fields 
we have to analyze the state of the matter field in the absence of matter, i.e., in 
the state of lowest energy. In an intuitive wave theory based on equations (4) 
and (6), the state of lowest energy is given when all electron states of negative 
energy are occupied and all states of positive energy are empty. In the presence 
of only a magnetic field, the stationary states of an electron can be divided into 
those of negative and positive energy. Hence the state of the lowest energy of 
the matter field can be derived in the same way as for a field-free space. 

The situation is different in an electric field. In this case the potential 
energy grows linearly in one space direction, so all energy values from — oo to 
+oo are possible. The eigenfunctions associated with different eigenvalues are 
transformed into each other by a spatial translation. Hence a classification of 
energy values into positive and negative is not unique. 

This difficulty is physcially related to the fact that in an electric field, pairs 
of positrons and electrons are created. The exact analysis of this problem was 
performed by Sauter 8 . 

8 F. Sauter, ZS. f. Phys. 69, 742, 1931 
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In Fig. 1, the potential energy V(x) and the lines V(x)+mc 2 and V(x) — mc 2 
are plotted against the coordinate (the electric field is parallel x-axis). The 
calculations of Sauter show that the eigcnfunction associated to the eigenvalue 
Eq, for example, is large only in the regions I and II. In the region II, they 
decrease exponentially. Therefore, a wave function that begins being large in 
region I decreases slowly in region III where the transmission coefficient through 
region II (which plays the role of a Gamow-wall) calculated by Sauter has the 

order of magnitude e~ Rc i £ w . If we define \<£k\ — ^^f- as the critical field 

strength, we can also write e e 7r . As long as |(£| <C |(£fc|, pair creation is 
so rare that it can be practically ignored. Then it must be possible to find 
solutions of the Dirac equation playing the role of cigenfunctions, which are 

large in region I but stay small of the order of e e T in region III. Conversely, 
it must be possible to find solutions which are large in region III and small in 
region I. After that we can characterize the state of the lowest energy by all 
electron states being occupied whose eigenfunctions are large only in region III, 
while the others are unoccupied. The energy density at x is calculated from the 
differences of the electron energies with Eq [compare Eq. (31)]. By switching off 
the electric field adiabatically, the so characterized state of the system goes over 
into the state of the field-free space, in which only the negativ-energy electron 
states are occupied. 

Our calculations follow those of Sauter. If an external magnetic field 2$ and 
an electric field € are present, both pointing in x-direction, the Dirac equation 
reads: 

lih d 6 d . s d 

\c dt e 1 1 1 Bx dy 

+ ^(ihj-- f.|SB|^-/?mcj V = 0. (19) 

The movement in y- and z-directions can be separated from that in x-dircction 
by the ansatz: 

I 

We introduce a new operator K by 
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K = +« s iJi A -(-■«,(— P( _ £.[S8|^_^- mc . (21) 



It satisfies: 



This equation can be understood as a wave equation determing eigenfunctions 
u n (y). If we set 

we find (22) which becomes essentially the Schrocdinger equation of the oszil- 
lator, 



4'JttX 



(H n (rf) is the n-th Hermite polynom). The eigenvalues are: 

em\n 



(n = 0, 1, 2, ■ • •) 

The operator K anticommutates with ot\ in the wave function (19), which can 
also be written as 

We can perform a canonical transformation in \, so that cr x becomes diagonal, 
and K and a\ read: 



m? c l + ■ 



{an + l + ffj 



| | 
— i of 



1 
— 1 



(26) 



Here the two matrices refer to another index independent of the spin orientation 
(they refer to the "g" -coordinate). One can treat a x as if it was an ordinary 
number (a x = ±1) and obtains with the convenient shorthand notation 



* = 



the equations 



(^ + i£)g + kf = o, 
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where / and g are the two components of the function x (associated with the 
"f/'-index). The equations (28) are formally identical with the Equations (12) 
from Sauter"s paper (op. cit. [HJ. They differ from those, however, in the mean- 
ing of k and that the system (28) actually has to be written down twice; once 
for a x = +1 and once for a x — —I. Sauter obtains two sets of solutions: 



* 2 . * 2 1 



1 « / ■ - * 2 > • * 2 1 



(29) 



The integrals are to be taken along a path coming from +oo circulating around 
both singular points +i/2 and — i/2 in a positive manner, and returning to +oo. 

Since we ignore pair creation in the calculation we consider as eigenfunctions 
only those parts of the functions / and g which vanish in one half of the space. 
Thus we set 

_ | U fiir £ > 0, 10 fur £ > 0, 

The new functions etc. do not correspond exactly to the stationary states 
but represent wave packets for which it is unlikly to diffuse into the initially 
empty area. For the construction of the density matrix we consider the states 
f\, g\, ... as being occupied, the states ff, g\, ... as unoccupied. We have 
doubled the number of states by the process (30). Thus, if we take all /-J, g{, 
/2, g\ as occupied and /-p, gf, /f , g\ as unoccupied, we shall obtain the density 
matrix twice. 

To find the energy density of the vacuum according to the method described 
in the first section we would have to calculate the density matrix associated with 
a finite distance of the two positions r' and r". From this we have to subtract the 
singular matrix S. Then obtain the energy density by taking the limit r' = r". 
For the following calculation it is more convenient to set from the beginning 
r' = r", but to restrict the summation over the stationary states to those of 
finite energy only. Aquivalently, we could make the summation convergent by 
an auxiliary factor e ~ conBt l E -( mc ) ], If the constant number in the exponent 
goes to zero, some of the terms in the energy density become singular, and they 
will be compensated by corresponding elements in the matrix S. The remaining 
regular elements yield the desired result. 

Bevor writing down the density matrix, we have to normalize the eigenfunc- 
tions. We can imagine the space of the eigenfunctions in x- and z-direction 
being confined to a large length L (the above eigenfunctions u n (y) are already 
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normalized). Then we have from the z-direction a normalization factor 1/y/L; 
from the x-direction, the asymptotic behavior of the Sautcr cigcnfunctions [com- 

fc 2 7r 

pare Eq. (22)] yields the factor 2-^=e~^. The summation over all states is to 
be done over all momenta which have the form 



Pz = -rm + const 
L 

and over all energies of the form E = jjj^m + const. These two sums can 
be converted into integrals. When neglecting the factor the differential is 

^r^Wc- -^ we calculate the density at position x , the energy of a state is the 
difference E — e\£\x - We therefore obtain for the energy density corresponding 
to the matrix Rs (compare Section 1) [for the meaning of a compare (33)] 



flP + WP-|/?P-|»f|'l,-(f-i) m 



which due to (23a) goes over into 



5}* S if) V -w ) e 2 



a 

— 1 



From this expression we have to derive the limit a— >0. We introduce the fol- 
lowing convenient abbreviations: 

These are dimensionsless and represent the ratio of the field strengths with the 
critical field strength \£k\, i.e., with the "l/137th fraction of the field strength 
at the radius of the electron" . With the equation (29) one finally finds 



-co 

< ■ 



H' 



— — log - 



fi ^&m*&4 v* 1 « ^ (34) 



where h(s\S2) is 
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The integration over £ yields: 



oe -f- 1 



p =-^S' o S o)a - fc - me *(x/\f ds 4 



'im^* ( Si + t)( S2+ I) 



ds a e 



(85) 



The first of the two s-integrals, i.e., the one over s\, can be performed by de- 
forming the integration path in such a manner that only a loop around the pole 
si = S2 — a is left. Replacing S2 by s = S2 — a in the result, one obtains: 



03 + 1 < a 



where 



\2 



■(4-7 



!?(<-«) + a(* + (37) 

The integration is to be carried out along a path shown in Fig. 2. The con- 
tour comes from +oo passes between the four poles, and returns back to +oo. 
Alternatively, one may integrate along the imaginary axis from +oo to — oo. 
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The main contribution to the integral comes from the regions between the 
poles. One may expand the logarithm in the exponent in power series of a, and 
obtains 



'.-(HP) 



For the following calculations, it is assumed that the electric field is small against 
the critical field \<&k\, i- e - «< 1 and thus 

fc' = Ip + b(2n+ 1 + <r,)]> 1. 

,2 

Then the expression — 4 ^" 1 in the exponent must be retained; the higher terms 
in the exponent may be considered as small, while we shall eventually take a— >0. 
They can therefore be included by an expansion. Thus, one obtains an expres- 
sion for f(k) of the form 



-i/2 



s 3 t 1 



(1 + 4s 2 )-^ + Bl? + Ck*. ..]. (39) 



Before performing the integral, it is convenient to calculate the summation 
over n and a; it is carried out according to the scheme 




It turns out that the higher terms of Euler summation formula contribute noth- 
ing to the end result. Finally, setting § = £, one gets for e — > (where 7 denotes 
the Euler constant 7 = 1.781). 
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(1 + 1nhf m ~i 



(The coefficients c m are to be determined further below). 

The parts corresponding to the singular matrix S still have to be subtracted 
from the result. One obtains easily the field- independent part of the singular 
energy density, which is to subtract by repeating the above calculation with 
plane waves: 



It is more difficult to calculate the field-dependent part of S. Since equation 
(13) contains only squares of field strengths in a and b, this would also for Us- 
Further, the constant C in equation (13) is chosen so that for constant fields, 
polarisation proportional to the field does not occur. From this it follows that 
the terms which are quadratic in the field are to be subtracted, and only the 
higher terms remain. This implies, anticipating an expansion in powers of b for 
b <C 1, that the first line of the right-hand equation (41) is to subtract on the 
whole. We proved this result by assuming that the points r/ and x" of the density 
matrix are different in the x direction [r' — r" = (a;, 0, 0)], and by setting at the 
end a = e = 0. After subtracting the terms which come from the matrix S, only 
the above discussed part of U remains. For the electric field, this calculation 
turned out to be very difficult. Starting with the total energy density, which is 
composed of the usual Maxwell energy density g^r(C£ 2 + 2J 2 ) and Dirac energy 
density U — Us, we go on with the Lagrangian using the relation (3) 
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One obtains 



(X)^116~B- 16 J -(l-21o g [l + (2^ + l)q) 



log [1 + (2 n'+l) 6] 



+ 

+ s 



i 1 +Znb 

2 re mi 



b 7 '' 

■ 2(l+2»6) Iog(l+2nfc) 
a i 

-log[l + (2 «' + !)&]) 



2m— 1 



(1 + 2m i») 2 



For small magnetic fields, the power series expansion of b may be found applying 
to equation (42) once more the Euler summation formula. Furthermore, since 
£ depends only on both invariants a 2 — b 2 and a 2 b 2 , from which it follows that 
£(a, b): £(a, 0) = £(0, ia), the missing coefficients c m , whose direct calculation 
would be very involved, may be determined indirectly from this relation. By 
calculating and C3 we have checked that the direct calculation of c m yields 
the same result; though we didn't find the general proof for it. In this manner, 
one obtains initially for small fields (a <C 1, b <C 1): 



;2 _ SB 2 ) + 16 7i>~ m c 2 jJM 



(a 2 — fc' J ) s + 7( 



+ 



13 (a &) 2 (a 3 



180 



The opposite limit (a< 1, &> 1) yields 

£ ^ I (<£ 2 - 93 s ) + 16 ri> m c* (^J {&» [1 log & — 0,191 j 



log 6 



+ T [log & - 0,145] + ^- + 0,202 - - [log 6 + 0,116] 

' + & [r2 + ^ + -] + -} 

We tried to derive an integral representation of £ in order to better see the gen- 
eral behavior of £ for an arbitrary field. This is possible when using the usual 
integral representation of the zeta function. One obtains 
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+ 





— ia brf - n ; , ; ' ' , r-p + 1+ ~ (i> 2 — « s ) ■ (45) 

' Cos (b + « a) »j — Cos (b — « a) tj 3 v J 

From the last form it is easy to see that £ only depends on both invariants 
l£ 2 — *8 2 and (£23) 2 . The cos terms allow an expansion in squares of the argu- 
ment (6 + ia) 2 — b 2 — a 2 + 2i(ab) and (b + ia) 2 = b 2 — a 2 — 2i(ab). Since the 
overall result is real, it may be represented as power series in b 2 — a 2 and (a6) 2 , 

which may be replaced by ^gji and jgjjx respectively. ^|<£fc| 2 = ^^f-^ • Thus 
the Lagrangian for arbitrarily oriented fields reads: 
2 = J p - 93 s ) 

- «-*!j£ ijfpB) J fe| — 

77 1 cos — 23 2 + 2 t (@ 23)) — konj. 

+ |«* | 2 4 £(8' -<£')}■ (45a) 

The first expansion term of equation (43) agrees with the results of Euler and 
Kockel (op. cit. EJ). 

The question of convergence of this power series expansion deserves closer 
investigation. The integral (45) converges for every value of b when a = 0. 
But when a ^ 0, the integral looses its meaning for rj = 221, • • • , where cot ar\ 
becomes infinite. According to this, the power series expansion in a, from which 
we started, can only be semiconvergent. We can give a definite meaning to the 
integral (45) by selecting an integration path that avoids the singular values 
— , — . But then the integral (45) acquires additional imaginary terms that can't 
be interpreted physically. Their meaning becomes clear if we estimate their size. 
The integral (45) around the pole r\=\ has the value — — • 4a 2 mc 2 ( I | £ ) 3 • e~ « 
(for 6 = 0). This is the order of the terms which are associated with the pair 
creation in the electric field. Thus, the integral is to be interpreted similar to 
the integration of the resonance denominator in perturbation theory. One may 
assume that a damping term corresponding to the frequency of the resonance 
process is responsible for the convergence. Further, one may assume that the 
result, which is obtained by avoiding the singular values, is correct except for 
the terms whose order corresponds to the frequency of the resonance process. 

The deviations from the Maxwell theory remain very small according to 
(43) and (44) as long as € and <8 are small against the electric field present at 

. 2 

a distance of x/137 • -^-r from the center of the electron. But also in the case 
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that the magnetic field exceeds this value, the additional terms to the Maxwell 
equations remain small (proportional to the original terms). They have the 
relative order j- as long as log b is near unity. Thus, the deviations from the 
usual Coulomb force between the protons (for example) come from the terms 
(43) and (44), and remain always very small. However, in this estimate, one 
must take into account that for a Coulomb field, the neglected terms containing 
the derivatives of the field strengths are more important than those contained 
in equations (43) and (44). 

3 Importance of the Result in the Quantum The- 
ory of Wave Fields 

The results derived in the last section may not be transfered immediately to 
the quantum theory of wave fields. In fact, one may show that the state of 
matter in a homogeneous field in the quantum theory of wave fields is not 
described by the above equations. Even if we start with the state of matter 
discussed in the last section as an "undisturbed state" , there are matrix elements 
of the perturbation theory which belong to the simultaneous formation of a light 
quantum and a pair. Even, if the energy does not suffice for a real formation of 
these particles, such matrix elements give rise to a perturbative energy of second 
order. This energy comes about through the virtual possibility of formation and 
disappearance of the light quantum and a pair, and the calculation provides a 
divergent result for them. The appearance of these perturbation terms may be 
made plausible by observing that the circular trajectories in a magnetic field are 
not really stationary states, but that the electrons in these states can radiate. 
In the classical theory of wave fields, this radiation can be ignored - this being 
crucial for the physical content of the calculation of the last section. Indeed, in 
the final solution the charge density and current density of matter vanish, so 
that radiation does not occur. In the quantum theory of wave fields on the other 
hand, this radiation does appear in form of the divergent perturbative energy 
of second order. 

One may also understand that a perturbative energy of exactly the same 
type occurs also in the field free vacuum ("self-energy of the vacuum"). Such 
self-energies are found always when calculating the contribution of second order 
or higher to the energy, which arises from virtual transition to another state and 
return to the initial state. These self-energies have so far always been ignored. 
For instance, the cross section of the Compton scattering may be calculated by 
performing a perturbation calculation to second order. Considering the terms 
of the forth order, one would get a contribution of the described type, which 
does not yield a convergent result. The calculation of the light-light scattering 
is carried out to the forth order (the lowest order that yields a contribution to 
the relevant process). The contributions of the sixth order would already di- 
verge. The previous success of these calculations, for instance the Klein-Nishina 
formula, seems to show that omitting the divergent contributions of the higher 
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order leads to the correct results. If this is the case, then it follows from the 
above discussion that the results of Section 2 can be transferred to the quantum 
theory of wave fields. This is also physically plausible, since the occurrence of 
the above radiation terms would remain incomprehensible according to the cor- 
respondence principle. Each single term in the expansion of the energy density 
in powers of € and 58 can concretely be associated with a scattering process, 
whose cross section is determined by it. The terms of forth order, for instance, 
stand for the usual light-light scattering, the terms of the sixth order determine 
the cross section of the process whereby three light quanta are scattered on 
each other etc. Irrespective of the question whether it is physically acceptable 
to neglect higher order terms, each expansion term in the result of the last 
section agree with a direct calculation of the corresponding scattering process 
in the quantum theory of wave fields if the perturbation calculation is only 
performed to the lowest order that yields a contribution to the corresponding 
process. In both calculations, the contributions of the terms are neglected which 
correspond to the formation and disappearance of the light quantum and a pair. 
[The agreement of the terms of forth order with the terms obtained by the di- 
rect calculation of light-light scattering (compare |3Jl is therefore a test for the 
correctness of the calculation.] For this reason, it is not ruled out that even 
the results for |*8| ^ |£fe| can be applied to experience. This is, however, cer- 
tainly not true for |€| > |€fc|; The pair creation destroys the basis for the above 
calculations when the electric fields are large. 

4 The Physical Consequences of the Result 

The results derived in the second section are very similar to Born's 9 derivation 
of Maxwell equations. Also Born obtains a more complicated function of the 
two invariants, € 2 — Q3 2 and (i£<8) 2 in addition to the classical Lagrangian 
2=2 _ sg2 incidentally, because of the actual value of §-, this function agrees 
with (43) in the order of magnitude of the lowest expansion terms (compare |3J). 
Though it is also important to emphasize the differences of both results. Born 
has chosen the modified Maxwell equation as the starting point of the theory, 
whereas this change in the Dirac theory is a really indirect consequence of the 
virtual possibility of the pair creation. In addition, Dirac's theory predicts 
also terms, which contain higher powers of the field strengths (compare ^ 121 
). Thus, especially the question of the self-energy of the electrons cannot be 
decided only with the help of these changes. The result of Born's theory imply 
that the changes of the Maxwell equations calculated here suffice to remove 
difficulties with the infinite self energy is an important indication of the further 
development of the theory. 

In this context, we must also ask the question whether the results about 
light-light scattering etc. derived from the Dirac theory may be considered to 
be final or whether it is expected that a future theory will lead to another 

9 M. Born, Proc. Roy. Soc. London (A) 143 , 410, 1933; M. Born and L. Infeld, Proc. 
Roy. Soc. London (A) 144 , 425, 1934; 147, 522, 1934; 150, 141, 1935. 
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results. Without doubt, the theory of the positron and the present quantum 
electrodynamics must be considered as temporary. Especially, the rules of find- 
ing the S matrix seem to be arbitrary in this theory (the inhomogeneity of 
Dirac's equation). The only reason for the matrix S introduced in 0] is the 
relative mathematical simplicity (together with some postulates on the formu- 
lation of the conservation laws). Thus, a deviation of the future theory from 
the present theory may well be possible. Since such variations may have a big 
influence on the Maxwell equations, the present theory is only reliable in its 
order of magnitude and the qualitative form of these changes. Up to now, it is 
hardly possible to make definite statements as to the final form of the Maxwell 
equations in future theory, since it is essential to consider all the processes of 
particles with very high energy (for instance, the occurrence of "showers"). 
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